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Transonic solutions of isothermal galactic winds in a cold 
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ABSTRACT 

We study fundamental properties of steady, spherically symmetric, isothermal galac- 
tic outflow in appropriate gravitational potential models. We aim at constructing a 
universal scale free theory not only for galactic winds, but also for winds from clus- 
ters/groups of galaxies. In particular, we consider effects of mass-density distribution 
on the formation of transonic galactic outflows under several models of the density 
distribution profile predicted by cosmological simulations of structure formation based 
on the cold dark matter (CDM) scenario. 

In this study, we have clarified that there exists two types of transonic solutions: 
outflows from the central region and from distant region with a finite radius, depending 
upon the density distribution of the system. The system with sufficiently steep density 
gradient at the center is allowed to have the transonic outflows from the center. The 
resultant criterion intriguingly indicates that the density gradient at the center must 
be steeper than that of the prediction of conventional CDM model including Navarro, 
Frenk & White (1997) and Moore et al. (1999). This result suggests that an additional 
steeper density distribution originated by baryonic systems such as the stellar compo- 
nent and/or the central massive black hole is required to realize transonic outflow from 
the central region. On the other hand, we predict the outflow, which is started at the 
outskirts of the galactic center and is slowly-accelerated without any drastic energy 
injection like starburst events. These transonic outflows may contribute secularly to 
the metal enrichment of the intergalactic medium. 

Key words: galaxies: evolution - galaxies: intergalactic medium - galaxies: ISM - 
cosmology: dark matter - ISM: jets and outflows. 



1 INTRODUCTION 

Galactic winds are widely thought to be an essential ingre- 
dient in the evolution of galaxies. It plays a key role by 
which energy and heavy elements are recycled in galaxies 
and are deposited into the intergalactic space. In the mod- 
ern paradigm of the galaxy formation based on the cold dark 
matter (CDM) hypothesis, it deduced that galaxies formed 
hierarchically in a bottom-up fashion, where a larger sys- 
tem results from the assembly of smaller dark matter halos. 
Baryonic gas falls into the gravitational potential well of 
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dark matter halos, and condenses rapidly as a result of the 
radiative cooling for atoms or molecules. The galactic winds 
seem to be inevitably influenced by the dark matter gravi- 
tational filed. We here study the fundamental nature of the 
galactic winds in the dark matter gravitational potential. 

The existence of heavy elements in the intra-cluster 
medium of galaxy clusters and the low-density intergalac- 
tic medium at z ~ 3 is the clear evidence of galactic out- 
flows (Songaila 1997; Ellison et al. 2000). Large number of 
large-scale outflows from star-forming galaxies has been ob- 
served in local galaxies (Lehnert & Heckman 1996; Martin 
1999). In the spectra of high- z galaxies, there are evidences 
for large-scale outflows. For instance, Lyo emission with red 
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asymmetric or P Cygni-type profile is commonly seen in 
z > 5 Lya emitters (Dey et al. 1998; Ellis et al. 2001; Daw- 
son et al. 2002; Ajiki et al. 2002). In addition, extended Lya 
nebulae, so called, Lya blobs at z > 3 have been possibly 
interpreted as spatially resolved large-scale winds driven by 
starbursts (Francis et al. 2001; Ohyama & Taniguchi 2004; 
Mori & Umemura 2006). Lyman break galaxies at z > 3 
also show signs of outflowing gas with velocity of several 
hundred km s" 1 . Furthermore, Weiner et al. (2009) and Ru- 
bin et al. (2010) recently reported that galactic outflows are 
ubiquitous in the optically selected star-forming galaxies at 
medium redshift z ~ 1 — 2 using stacked spectra from the 
DEEP2 survey. 

In the theoretical point of view, so far, the large num- 
ber of studies about the galactic winds has been produced 
by numerical simulations. The recent developments in com- 
puter technologies and numerical methods have made it pos- 
sible to simulate the multidimensional dynamical or chemo- 
dynamical evolution of galaxies including the effect of star 
formation and supernovae (SNe) feedback. For example, the 
simulation of the galactic winds from elliptical galaxies are 
motivated by the observation that elliptical galaxies have 
very little interstellar medium. While the numerical simula- 
tion has become a powerful tool to explore the formation and 
evolution of galaxies, a simple analytical approach plays an 
indispensable role to reveal the essential nature of the ac- 
celeration process of winds and to interpret observational 
data. 

Fundamental nature of thermally driven spherically 
symmetric transonic outflow was examined by Parker (1958) 
in relation to the solar wind, and then it has been dis- 
cussed for more than 50 years. This model was based on 
a hydrodynamic description of the sun's atmosphere, and 
proposed that the solar wind is a smooth, spherically sym- 
metric, time-steady transonic outflow of hot gas. We must 
note that the transonic solution is very special and plausible 
one, because it is the entropy maximum solution connecting 
the sun and infinity. Thus, such transonic solution is widely 
accepted as the most plausible universal solution of the out- 
flows from any object. Parker's solution shows that the flow 
inside the transonic point resembles the hydrostatic equilib- 
rium closely, whereas the flow tends to free-expansion out- 
side the transonic point. Then, Holzer and Axford (1970) 
provided the generalization of Parker's approach to galac- 
tic outflow (see also Burke 1968; Chevalier & Clegg 1985). 
Wang (1995a, 1995b) studied for a radial, steady solution 
of galactic wind including the effect of the galactic gravita- 
tional potential and the efficiency of radiative cooling. As- 
suming a single power-law distribution of mass density, they 
focused on the outflow that is supersonic everywhere from 
less massive galaxies. Recently, Everett & Murray (2007) 
investigated the Parker-type acceleration of galactic wind 
under realistic but complicated situation with gravity from 
a point mass, cloud/ISM drag force, adiabatic cooling pro- 
cess, and photo ionization/absorption process. They con- 
cluded that the outflow accelerated insufficiently under adi- 
abatic cooling. Moreover, Sharma & Nath (2012) emphasis 
that the importance of extra energy source such as photo 
ionization/absorption process by AGN radiation or SNe en- 
ergy/momentum injection in order to attain sufficient accel- 
eration. 

However, recent cosmological JV-body simulations based 



on collisionless CDM scenario no longer indicated such a 
single power-law distribution but always predicted a dou- 
ble power-law density distribution. Navarro, Frenk & White 
(1997) pointed out that their structure can be approximated 
as Pt>m(v) oc r~ x (r + ro)~ 2 , where r is the radius from the 
center of the galaxy and ro is the scale radius at which 
the density profile agrees with the isothermal profile (i.e., 
Pdm(^) oc r~ 2 ). This is related to the formation epoch of the 
CDM halo (see also Navarro, Frenk & White 1997). Then, 
Fukushige & Makino (1997) and Moore et al. (1999) used 
high-resolution simulations and showed that the CDM halos 
have a steeper central cusp than quoted above. The result- 
ing structure of the CDM halos depends on the number of 
particles used in the simulation and is still an open question 
(see Navarro et al. 2010; Ishiyama et al. 2011). In contrast, 
recent observations of nearby dwarf galaxies and low surface 
brightness galaxies have revealed that the density profile of 
the dark matter halo is constant at the center of such galax- 
ies. For instance, Burkert (1995) proposes that the density 
profile, pdm(»*) oc (r + ro) (r 2 + r§) , nicely reproduces 
the rotation curves of nearby dwarf galaxies and the central 
density is correlated with ro through a simple scaling rela- 
tion (e.g., Moore 1994; Burkert 1995; de Blok et al. 2001; 
Swaters et al. 2003; Gentile et al. 2004; Spekkens, Giovanelli 
& Haynes 2005). This is well-known as an unsolved problem 
in the CDM scenario, and so-called the 'Core-Cusp prob- 
lem'. 

Under these inconclusive circumstances for determining 
CDM distribution, no model considers facing the transonic 
galactic winds in a realistic constellation of the CDM distri- 
bution of galaxies. These situations motivate us to explore 
the series of the solution for the transonic galactic winds in 
the appropriate CDM halo models. In this paper, we focus 
on the analytical steady-state solution for a spherically sym- 
metric isothermal outflow in a various CDM halo models. 
We aim at clarifying influence from the dark matter gravi- 
tational field on the galactic winds. We are interested espe- 
cially in realizability of the transonic winds. This research 
is conducted to be a scale free universal wind theory and 
does aim not only at winds from galaxies but also at winds 
from groups/clusters of galaxies because CDM halo density 
distribution model seems to be rather universal for these 
hierarchies of objects (Navarro, Frenk & White 1997). For 
analytical convenience, we assume isothermal flow without 
any mass injection along the flow, because gas temperature 
in actual galaxies seem to be constant (Li et al. 2011) and 
the locus of the resultant transonic point in this study is far 
distant from the star forming region (see Fig. [2] of this pa- 
per). The adequacy of these assumptions is discussed later. 

The structure of this paper is as follows. In Section 2, we 
describe the basic equations for the analysis of the transonic 
galactic winds. In Section 3, we show the critical condition 
for the existence of the transonic winds and their series of 
solutions. Finally, we discuss the results in Section 4. 



2 TRANSONIC FLOW IN CDM HALO 

2.1 Basic equations 

We consider steady, spherically symmetric, isothermal gas 
outflows from galaxies. The basic equations for this problem 
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are mass and momentum conservation laws: 



4nr pv — M — const, 



and 

dv 
dr 



c s dp d(j> 
p dr dr ' 



(1) 



(2) 



where p,v,c B ,M, and <f> are the density, the velocity, the 
sound speed of the gas, the mass-loss late, and the gravita- 
tional potential of a galaxy, respectively. We introduce the 
nondimensional distance x = r/ro from the galactic center 
where ro is the scale radius of the galaxy. Substituting p 
from equation (fl]) into equation {5), we obtain 



c 2 dM _ (2c?/x) - (dcj)/dx) 



M dx 



M 2 



(3) 



where A4 = v/c s is the Mach number. The numerator of 
the right-hand side denotes a change of the effective cross 
section of a Laval nozzle. Thus, we can see that the grav- 
itational force plays a role of choking of the cross section 
(hereafter, the gravitational choking). Integrating equation 
0, we obtain the relation, 



M 2 -InM 2 =41nx-^| + C, 

'■■ s 



(4) 



where C is the integration constant. A series of wind solu- 
tions is obtained from equation Q as M 2 (x) parameterized 
byC. 

The feature of solutions of equation Q is similar to 
Parker's solution of solar wind (Parker 1958). The denomi- 
nator in the right-hand side of equation ^ is negative for 
subsonic (M < 1), and positive for supersonic (M > 1). 
If the gravitational choking is effective (2c 2 /x < d(j>/dx), 
the numerator is negative, while if the gravitational choking 
is not effective (2c 2 jx > d<j>/dx), the numerator is posi- 
tive. It is equivalent to the cross section of a Laval nozzle 
that is decreasing (increasing) for negative (positive) nu- 
merator. Subsonic flow accelerates in the region where the 
numerator is negative and supersonic flow accelerates in the 
region where the numerator is positive. The singular point 
at which both the denominator and the numerator equal 
to zero (M — 1, 2c 2 jx — d(f>/dx), and sign of the numera- 
tor changes from negative to positive is called the transonic 
point(or the critical point, see Bondi 1952; Parker 1958). 
Only the transonic flow can continuously accelerate from 
subsonic to supersonic by passing through the transonic 
point. 



2.2 Density distributions 

The NFW model that is widely accepted density profile of 
CDM halos was empirically derived from cosmological N- 
body simulation (Navarro, Frenk & White 1997). Fukushige 
& Makino (1997) and Moore et al. (1999) derived steeper 
distribution from JV-body simulation with higher resolution. 
The resultant empirical profiles of density distribution is ap- 
proximately fitted by 



PDM(r;a) 



poro 



r a (r + r ) 3_ 



(5) 



where po is the scale density. In this equation, a = 1 is 
the NFW model and a — 1.5 corresponds to the model 



advocated in Fukushige & Makino (1997) and Moore et al. 
(1999). On the other hand, observed density profile derived 
by Moore (1994) and Burkert (1995) is represented by a = 0. 
The plausible value of the index a remains an open question. 
Thus, we treat a as a parameter and study the variation of 
solutions depending on a. By choosing ro as the unit of 
length, equation ((5| is rewritten as 



p D m(x;a) 



Po 



*(x + l) 3 



(6) 



where x = r/ro. In the limit of x — > 0, pom(x;a) oc x~ a 
and pT>u(x\a) oc x~ 3 for x — > oo. Note that the functional 
form of the gravitational potential near the center crucially 
depends on the index a. 



3 SERIES OF GALACTIC WIND SOLUTIONS 

3.1 Condition for transonic flow from the center 

The acceleration of outflows obviously depends on the shape 
of the gravitational potential from the discussion in section 
2.1, and the existence of transonic solutions starting from 
the center, therefore, relates closely to the power-law index 
a in equation {5). Here, we consider the critical condition 
upon q for the transonic solution from the center. 

The total mass M(x) within x under the density distri- 
bution ((6| can be described analytically using Gauss's hy- 
pergeometric function (see appendix A) as 



M(x;a) 



47rp r 3 _ 



2F1 [3 - 



,3-a,4- 



(7) 



if ^ a < 3. Substituting d<j)(x\a)/dx = GM(x;a)/(r x 2 ) 
into the numerator of the right-hand side of equation Q, 
we obtain 



2c: 



d(f>(x;a) _ 2c 2 
dx x 



i-E^f- 1 )" 3 



(8) 



where {A n } are positive constants (see appendix A). In the 
limit of x — > 0, equation ([8]) approaches 



2ci 



d<f)(x; a) 
dx 



2c: 



1 - A x z 



(9) 



If the choking by the gravitational force is effective in the 
innermost region, the numerator is negative. In the case of 
^ a < 2, the numerator is positive because l—Aox 2 ~ a — > 1 
as x —¥ 0. In the case of a = 2, 1 — Aox 2 ~ a — 1 — Aq. 
Therefore the numerator is negative if Aq > 1. In the case of 
2 < a < 3, the numerator is negative because 1 — Aox 2 ~ a — > 
—00 as x — > 0. Finally, we conclude that the condition for 
realizing transonic flow from the central region is a ^ 2. 

We note that this is rather steep density profile in the 
innermost region. Widely accepted values of a for the den- 
sity distribution profile of CDM halos are a = 1 and 1.5, but 
they do not fulfill the critical condition a ^ 2. The galactic 
wind solutions for the case of ^ a < 2 (Case 1), and the 
case of 2 Sj a Sj 3 (Case 2) are shown in the following sec- 
tions. Here, we focus only on the cases a = 0, 1 and 1.5 in 
Case 1, and a — 2 and 2.5 in Case 2 because of analytical 
simplicity. 
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3.2 Case 1: < a < 2 

3.2.1 The galaxy mass, the gravitational force, and the 
gravitational potential 

The total mass M(x; a) within x is given for specified values 
of a as follows, 



M(x; a — 0) = inporo 

M(x; a — 1) — 4nporo 

and 

M(x; a = 1.5) = 8nporo 



]n(x + 1) - 
In(a; + 1) - 



x(2 + 3x) 



2(x + l) 2 

X 


x + 1 


•) 



(10) 

(11) 



ln(>/S + \/iTT) - 



x + 1 



(12) 



Therefore the corresponding gravitational forces are given 
by 

#(I;Q = 0) =4WG 



dx 



ln(x + 1) 



d(b(x:a — 1) , 2,-, 1 

da; x A 

and 

d(f>(x; a = 1.5) 



x(2 + 3x)~ 




2(x + l) 2 _ 


> 


ln(a! + 1) 


X 

x + 1 



(13) 
(14) 



dx 



S%poro G 



ln(yi + Vx + 1) 



x + 1 



(15) 



respectively. By integrating equations (|13|l . (|14[) and (|15|) . 
corresponding gravitational potentials are given as 



(x; a = 0) = — Anporo G 



- hl(x +1) ; 

x y ' 2(x + l) 



<j>(x; a — 1) — —4-Kporo G— ln(x + 1), 
x 

and 

(j>(x;a = 1.5) = — 87rporo G 

1 

X 
respectively. 



(16) 

(17) 



[l + ln{V^+VxTT- \fxVxTT)], (18) 



3.2.2 Mechanism of acceleration 

From the numerator of right-hand side of equation ((3j), we 
define 



2c 2 A<j)(x\a) _ 2c\ 



N(x;K,a) 



(19) 



iV(x;J,,a)^l-J,x 2 - 2Fl[3 ~ a ' 3 ~ a;4 ~ Q; ~ xl , (20) 

3 — a 



and 



K 



2npor G 



Q 



(21) 




Figure 1. N(x; K, a) with a = 1: If K > K c , two points (A, B) 
at which N(x\K,a) vanishes appear. The gravitational choking 
is effective in x\ < x < xb (2:a,£b are the loci of A and B, 
respectively). This pattern does not topologically change in ^ 
a < 2. 



Substituting equations (|13[) , (|14|) or (JTSj , in sequence, into 
the numerator of the right-hand side of equation ([3]), we 
obtain 



N(x;K,a = 0) = 1 

N(x;K,a= 1) = 1 

and 

N{x;K,a = 1.5) = 1 



K 

X 

K 

X 



ln(x + 1) 



x(2 + 3x) 



2(x + l) 2 



ln(x + 1) - 



x + 1 



(22) 
(23) 



2K 



ln( Vx + \/xTT) 



x + 1 



(24) 



for each value of a. If the factor iV(x; K, a) is positive, sub- 
sonic flow decelerates and supersonic flow accelerates, while 
if N(x; K, a) is negative, subsonic flow accelerates and su- 
personic flow decelerates. Fig. [T] shows the behavior of the 
factor N(x; K, a) with a = 1 for example. We can find that 
there are three cases depending on K. A and B in Fig. [T] are 
the points at which N(x; K, a) vanishes. There is a criti- 
cal value K = K c such that minimum of N(x; K c ,a) is zero. 
The value of K c depends on a as follows: K c ~ 4.6 for a — 1, 
K c ~ 2.9 for a = 1.5, and K c ~ 8.2 for a = 0. If K < K c , 
N(x;K,a) is always positive. If K > -Kc, two points (A, B) 
at which A^(x; if, cv) vanishes appear. We name the loci of A 
and B as xa and xb, respectively, xa is the locus that sign 
of the factor A^(x; K, a) changes from positive to negative. 
xb is the locus that sign of the factor N(x;K,a) changes 
from negative to positive. Thus, the point (x, M 2 ) = (xb, 1) 
at which galactic wind passes xb with sound speed become 
the transonic point. 

3.2.3 Topology of the solution 

We substitute equations (|16|) . (|17|) . and (|18[) into equation 
([4]) successively and draw solutions in the plane (x, M 2 ) for 
each case of a. Fig. [5] shows family of solutions for a — 1 
for example. Fig. [2] is the so-called 'phase diagram' of the 
solutions of equation Q. We can see that the topological 
property of the phase diagram strongly depends on K as 
shown in Fig. Q] Fig.[2»,[2I>) and^c) are for K = 3.8(< 
K c ), K = K c ~ 4.6 and K = 7.5(> K c ), respectively. 
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Figure 2. Family of solutions for a = 1: The upper panel is 
for K = 3.8(< K c ), the middle panel and the lower panel are 
for K = K c ~ 4.6, and for K = 7.5(> K c ), respectively. In 
lower panel, sign of velocity gradient changes at x\ ~ 0.43 and 
xb ~ 12.7 and transonic solution which passes the transonic point 
(x,A4 2 ) ~ (12.7, 1) appears. The patterns of solution curves for 
^ a < 2 do not topologically change. 



If K ^ K c , subsonic flow decelerates and supersonic 
accelerates monotonically (Fig. HJa) and[2]D). Consequently, 
transonic solution does not appear in these cases. If K > K c , 
the subsonic flow accelerates and supersonic flow decelerates 
in x A < x < x B (Fig. 0:). The point at (x,M 2 ) — (x A , 1) 
is so-called 'O-point' from its topological property in Fig. 
[2fc). Similarly, the point at (x,M 2 ) — (xb, 1) is so-called 
'X-point'. This is the transonic point. In Fig. [2Jc), the sign 
of velocity gradient changes at x\ ~ 0.43 and xb — 12.7. 
The transonic point forms at (x, M 2 ) ~ (12.7, 1). Therefore 
a transonic solution which passes (x,A4 2 ) ~ (12.7, 1) exists. 

3.3 Case 2: 2 ^ a < 3 

3.3.1 The galaxy mass, the gravitational force, and the 
gravitational potential 

The total mass M(x; a) within x for a = 2 and 2.5 are given 
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Figure 3. N(x; K, a) with a = 2 (upper panel) and with a = 2.5 
(lower panel), respectively. Lower panel is for K = 0.50. The 
pattern of lower panel does not topologically change in 2 ^ a < 3. 



and 



M(x; a = 2.5) = 87rp r In (^/x + Vx + l) 



(26) 



respectively. The gravitational forces for each case, there- 
fore, are 



Ad>(x\ a — 2) , 2 _ 1 , . „. 

-Z±-L >- = 47rp ro G-5- In i + 1 , 

da; x z 

and 



2 51 1 

— - = 8TTp r 2 G— In {y/x + Vx + 1) 



dx 

respectively. By integrating equation (|2Tf) and 
itational potentials for each case are given as 

4>{x;a = 2) = — 47rporo G 



■ \n(x - 



In- 



x + 1 
and 
4>{x; a — 2.5) = — 8nporo G 

y/xy/x + 1 + In (y/x + \/x + 1) 

X 

respectively. 

3.3.2 Mechanism of acceleration 



(27) 

(28) 
the grav- 

(29) 



- 1 



(30) 



Substituting equations (I27[) or (1281) into the numerator of 
the right-hand side of ©, we obtain 



M(x; a = 2) = 47rp n) In(a: + 1), 



(25) 



N{x;K,a = 2) 



K 



ln(s + l), 



(31) 



6 Tsuchiya, Mori and Nitta 



and 



N(x;K,a = 2.5) 



2A 



In (y/x + Vx + T) 



(32) 



Fig. [3ja) shows N(x; K, a) for a — 2 and Fig. [3jb) shows 
a — 2.5 as a function of x. We can see that the function 
N(x; K, a) monotonically increases for a ^ 2: it is critically 
different from Case 1. In the Case 2, one X-point without 
O-point appears. In this case, we define B as shown in Fig. 
[3] In the case of a = 2, the value of N(x; K, a) is finite at 
x — and there is a critical value K = K c — 1. Therefore if 
K ^ 1, JV (x; A", a) does not become negative, while if A' > 1, 
N(x; K, a) becomes negative in x < xb (ib is the locus of 
B) in Fig. [3Ja). Thus, the point where (x,M 2 ) — (ib, 1) 
becomes the transonic point. In the case of a = 2.5, the 
value of the function is negative infinity in the limit of x — > 
0. Therefore a transonic point always exists at (x,M 2 ) = 
(z B ,l)mFig.[Ib). 



3.3.3 Topology of the solution 

By substituting equations (|29l) or (|30l) into equation @, we 
can draw phase diagram for each case. Figs. U and [S] are 
phase diagram for a — 2 and a — 2.5, respectively. Figs. 
[!», (b) and (c) are for K = 0.75(< K c ), K = 1(= K c ), 
K = 1.25(> A c ), and Fig. O is for A = 0.50, respectively. 
In the case of a = 2 (Fig. BJ, if A" > K c , xb at which 
the sign of the velocity gradient changes appears and there 
is a transonic solution which passes the transonic point at 
(x, M 2 ) = (xb, 1). In the case of a — 2.5 (Fig. [5]), there is 
also a transonic solution which passes the transonic point at 
(x,M 2 ) — (xb, 1). In this case, according to Section 3.3.2, 
we find such critical solution for any value of K. According 
to Figs. HJc) and \S\ the transonic wind can start from the 
center x = in Case 2. 



4 SUMMARY AND DISCUSSION 

The property of a stationary, spherically symmetric and 
isothermal galactic wind in a CDM halo is examined in this 
paper. Depending on the mass-density profile of the CDM 
halo, three types of outflow can exist: (i) supersonic or sub- 
sonic flows everywhere, and flows with (ii) an X-point or 
(iii) a pair of an O-point and an X-point. The condition for 
existence of a transonic solution is sensitive to the mass- 
density profile of the CDM halo. Especially, the transonic 
outflow from the galactic center is realized under necessary 
condition that the power-law index of the mass-density dis- 
tibution near the galactic center must be steeper than two 
(q ^ 2 in equation [S]) . 



4.1 The locus of the transonic point 

Our analysis shows that the locus of the transonic point is 
critically related to both the power-law index a of the mass- 
density distribution and the coefficient K (see equations [S] 
and !21|) . In Fig. [6] we summarize the locus xb of the tran- 
sonic point obtained by solving equation N(xb',K, a) = 
for ^ a < 3 and 1 ^C K ^ 19. Each solid-curve shows the 
locus xb as function of a for given value of K in the range 
from K = 1 (bottom) to 19 (top). Using equation Q and 



M< 



M< 



M< 




X 

Figure 4. Family of solutions for a = 2: The upper panel is 
for K = 0.75 < K c , the middle panel is for K = K c = 1, and 
the lower panel is for K = 1.25 > K c . In lower panel, sign of 
velocity gradient changes at xb — 0.54 and transonic solution 
which passes the transonic point (x,A4 2 ) ~ (0.54, 1) appears. 



M< 




x 

Figure 5. Family of solutions for a = 2, for K = 0.50: Sign 
of velocity gradient changes at xb — 0.81 and transonic solution 
which passes the transonic point (x, Ai 2 ) ~ (0.81, 1) appears. The 
patterns of solution curves for ^ a < 2 do not topologically 
change. 
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Figure 6. Locus of the outer X-point x c obtained by solving 
equation N(x c ; K, a) = for ^ a < 3 and 1 < K < 19. Each 
solid-curve corresponds K = 1, 3, 5, 7, 9, 11, 13, 15, 17 and 19 from 
bottom to top. Each dashed-curve corresponds K = /(o)/7 with 
7 = 1, 0.5 and 0.1 from bottom to top. 



Note that the virial radius is comparable to the locus of the 
transonic point for a = 1 in Fig. 6. Depending on the coef- 
ficient K and the virial mass of the galaxies, therefore, the 
wind is accelerated mainly in the region far outside the opti- 
cally visible scale of the galaxies in some cases. Moreover, we 
must note that the density structure of the transonic wind 
in the subsonic region x < xb is crucially similar to that 
of the hydrostatic equilibrium except in the vicinity of the 
transonic point. Thus, we may observationally confound the 
slowly accelerated wind structure having xsro ^> r v [ T with 
the hydrostatic one. 

There possibly exists such slowly accelerated outflows 
in the outskirts of galaxies with no drastic heating by such 
as starburst events but by quasi-stationary heating from the 
heat reservoir of shock heated CDM halo (see next subsec- 
tion). Such slowly accelerated galactic wind may play an 
important role of the metal enrichment of the intergalactic 
medium though it resembles closely to the hydrostatic gas 
in the optically observable region. 



Table 1. Numerical values of /(a) for typical a 



a 





0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


/(«) 


14.7 


8.86 


5.18 


2.87 


1.44 


0.567 






Cg = kBT/fj,m p , where fj, is the mean molecular weight, m p 
is the proton mass, fee is Boltzmann's constant, and T is the 
gas temperature, equation (|21|l becomes 



GM(r ;a)fim p 
xi = — ; — — j (a 



(33) 



2r k B T 

where M(ro; cv) is the enclosed mass within the radius r — ro 
for a given a, and f(a) is defined by 

3 



/(«) 



(34) 



2 Fi[3-a,3-a,4-a;-l]' 

The numerical values of f(a) for typical a are summarized 
in Table [T] Accordingly, the coefficient K is roughly inter- 
preted as the ratio of the gravitational potential energy of 
the galaxy to the thermal energy of the gas. 

Fig. [B] indicates the allowable range of the transonic 
point for given K. For instance, the transonic solution can 
exist only in the range 2 ^ a < 3 for K = 1 and ^ a < 
3 for K > 8.15. Especially, the transonic solution starting 
from the center x — can exist only for a ^ 2. This is 
consistent with the discussion in section 3.1. In addition, we 
can see that the locus of the transonic point is an increasing 
function of K for given a. For the NFW model (cv = 1), the 
transonic point locates at Xb — 4.22 for K — 5 and xb = 
21.6 for K — 10. Incidentally, the virial radius of the CDM 
halo is usually defined by r v i r = cro for the NFW model, 
where c is so-called the concentration parameter. Using the 
concentration-mass relation 



log 10 c = 0.971 - O.O941og lo (M/lO 12 /i" 1 M ), 



(35) 



derived by Maccio, Dutton & van den Bosch (2008), we ob- 
tain c as a function of the virial mass. This equation gives 
7 i$ c < 30 with respect to a reasonable range of the virial 
mass 7 < log 10 (M/M Q ) < 13 for the NFW model. In other 
words, the virial radius is in the range 7ro < r v j r < 30 ro. 



4.2 Availability of our simplified model 

4-2. 1 Isothermal approximation 

In the standard CDM picture of galaxy formation, gas falling 
into CDM halos is shock-heated approximately to the halo 
virial temperature, thus the CDM halo maintains quasi- 
hydrostatic equilibrium (Rees & Ostriker 1977; Binney 1977; 
Silk 1977; White & Rees 1978). We assume that the wind 
gas is in thermally equilibrium with such heated CDM halo. 
In this model, CDM halo plays a role of a heat reservoir to 
keep wind gas temperature nearly to be the halo virial tem- 
perature. In this case, the gas could be treated as isothermal 
and its temperature T is expected to be close to the virial 
temperature in a wide range up to the virial radius, 



Tvir = 7 



Gfj,m p M(r ; a) 
2k B r 



(36) 



where 7 is a fudge factor of the order of unity that should be 
determined by the efficiency of the shock heating. A typical 
value of T v i r is ~ several xlO 6 K for galaxies with mass ~ 
10 11 M and scale radius ~ 10 kpc. 

Indeed, in case of the Sombrero galaxy, the gas temper- 
ature seems to be isothermal in a considerably wide spatial 
range ^ 25 kpc from the galactic center (Li et al. 2011). We 
note that the temperature of the gas in the Sombrero galaxy 
(~ 0.6 keV ~ 7 x 10 6 K) is consistently recognized as the 
virial temperature (see also Mathews & Brighenti 2003). On 
the other hand, recent X-ray observations reveal a variety 
of the temperature profiles of the hot interstellar medium 
in early-type galaxies. Diehl & Statler (2008) categorized 
the observed temperature profiles into four major groups: 
isothermal (flat), positive gradient (outwardly rising), nega- 
tive gradient (falling) , and hybrid (falling at small radii and 
rising at larger radii) (see also Fukazawa et al. 2006). We 
have, however, a poor understanding the origin and evolu- 
tion of these temperature gradients. In addition, since X-ray 
emissivity in the outskirts of the galaxies is very low to be 
detected with current X-ray satellites, the maximum radius 
of the observed X-ray emission is usually smaller by a fac- 
tor of 10 than the virial radius of the CDM halo (Fukazawa 
et al. 2006). Thus, there is a significant ambiguity keeping 
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the temperature gradient in the outskirts of the galaxies. 
Under these complicated situations, we focused on the flat 
temperature profile in this paper. 



4-2.2 Is mass injection along the wind necessary? 

One may think that mass injection along the wind flow due 
to the SNe is essential and cannot be neglected for plausible 
galactic wind models. Actually, it may play intrinsic roles as 
material/energy source of the wind for central star forming 
region of galaxies that is in a locus of smaller radius than 
the transonic point. In addition, mass injection results in 
effective braking force on the wind like gravity. These make 
essential influences on the wind nature in the central region. 

However, in the resultant wind solution of our analy- 
sis for actual CDM halo models, the transonic point forms 
at far distant region from the galactic center (~ 10 times 
the scale radius of the CDM halo density distribution that 
is much larger than the radius of typical locus of the star 
forming region, see Fig. [2}. In such the distant region, we 
can neglect any effect from star formations, thus can neglect 
mass injection by SNe along the wind flow when we discuss 
the acceleration process of the wind. Moreover, when we ap- 
ply our result to the wind from groups/clusters of galaxies, 
the mass injection along the flow also can be neglected. Any 
effect from mass injection is meaningful only at the region 
in the vicinity of the starting point of the flow, and thus, it 
does not play any important role on the actual acceleration 
process of the wind if the transonic point forms sufficiently 
outside the scale radius as like our case. 

We can conclude that the mass injection along the wind 
is not necessarily important for galactic wind theory except 
the case that the transonic point is very close to the galactic 
center. This simple but analytically effective approximation 
we adopted in this paper may be allowed as the first step to 
explore the transonic nature of the galactic winds. 



Rees 1992; Efstathiou 1992). In this case, K is no longer the 
function of the one-parameter family of a quoted above. 



4.3 Implication of the critical condition a ^ 2 

We have revealed that transonic solution starting from the 
galactic center needs rather steep density distribution (see 
Sec. 3.3). This is a natural consequence that the solution 
must approach asymptotically to the Parker's solution in 
the limit of steep density gradient. It is striking that none of 
cosmological iV-body simulations based on the collisionless 
CDM predicts such a steep power-law index (a ^ 2) around 
the center of CDM halos. Therefore, the density enhance- 
ment of the baryon at the central region of the CDM halo 
may be essential for the transonic galactic wind from the 
galactic center. In other words, the gravitational potential 
induced by a stellar system and/or a central massive black 
hole plays a crucial role for the acceleration of the galactic 
wind from the galactic center. 

Though, we studied only the effect from the CDM halo 
as a source of gravitational potential in this paper, it is very 
interesting to examine the case including the gravitational 
potential induced by the stellar system and the central mas- 
sive black hole. We would then also need to consider the 
effects of the star formation and subsequent feedback pro- 
cess such as stellar winds and SNe heating from massive 
stars. These feedback processes will supply thermal energy 
into the gas and will accelerate the outflows. 

In a series of forthcoming studies, we plan to report 
the results taking into account the multiphase states of the 
gas with cool components in a galaxy, including the effect 
of the radiative cooling of the gas, energy input from stars 
and AGN, and the UV background radiation. In this case, 
the efficiency of the acceleration of the outflow may be quite 
different. 



4-2.3 Range of the coefficient K 

We here discuss the range of expected value of the coefficient 
K defined in equation l21l As discussed in Sec. 3, the value of 
K is critical to determine the topology of the solution curves 
in the M 2 — x phase diagrams (see Figs. 1-5). However, note 
that the value of K is still uncertain both observationally 
and numerically. 

Substituting T = T v i r into the equation (|33[) . we obtain 



K = /(a)/ 7 . 



(37) 



Eke, Navarro & Frenk (1998) adopted 7 = 2.3 and Kitayama 

6 Suto (1997) adopted 7 = 1.8 as their canonical value in 
the analysis of galaxy cluster number counts. For instance, 
assuming 7 = 1, NFW model (a — 1) yields K — 5.2 and 
Fukushige-Makino-Moore model (a = 1.5) yields K — 2.9. 
According to Table [l] a reasonable range is < K < 15 in 
this case. Each dashed-curve corresponds K — /(a)/7 with 

7 = 1, 0.5 and 0.1 from bottom to top in Fig. [5] 

In addition, the combination of the extra heating mech- 
anisms such as SNe and UV background radiation and the 
radiative cooling process could keep the temperature to be 
constant in many situations (see Dekel & Silk 1986; Yoshii 
& Arimoto 1987; Mori, Yoshii & Nomoto 1999; Babul & 
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APPENDIX A: GAUSS'S HYPERGEOMETRIC FUNCTION 

Hypergeometric function defines as follows. 

(ai)n ■ ■ ■ (flp)n Z n 



p F q [ax,--- ,a P ;bx,- ■■ ,b q ;z] = ^ 



(6l)n ■ ■ ■ Q>q)„ n! 



(a)„ — a(a + 1) ■ ■ • (a + n — 1) = 



T(a + n) 
T(a) ' 



(o)o = 1. 

Especially, in case of p = 2, g = 1, it is called Gauss's hypergeometric function. It is given by 

T(c) ^ r(o + n)T(b + n) 2 n 



2 Fi [a,6;c;z] 



r(o)r(6) 



E 



T(c + n) n! 



(Al) 

(A2) 
(A3) 

(A4) 



When using a directly, without substituting a number for it, total mass of galaxy within x, M(x) is given only by series 
expansion, using the Gauss's hypergeometric function. It is given by 



M{x) = 



47rp ro 3 



3-a 



x a 2F1 [3 — a, 3 — a; 4 — a; —a;] . 



Substituting equation (|A5|) in numerator in the right side of equation @, we obtain 

(2c s 2 /:r) - (d^/dx) = — - 47r o P ° r ° ^ '^ 2F1 [3 - a, 3 - a;4 - a; -x] , 
X A — a 



2c 2 s 4irpor 2 G T(4 - a) 
x 3- a r(3-a)r(3-a) 



E 



r(3 - a - n)r(3 - a - n) 1 
T(4 — a — n) n\ 



(-l)"i" +1 " a , 



'2c! 



1 - E A n (-l) n x n+2 - a 



2TTp rg 2 G T(4 - a) r(3 - a - n)r(3 - a - n) 

c|(3-a) r(3-a)r(3-a) r(4-a-n) ' 



(A5) 

(A6) 
(A7) 

(A8) 
(A9) 



